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Abstract 

Based on the recently developed light-front similarity renormalization group 
approach and the light-front heavy quark effective theory, we derive ana- 
lytically from QCD a heavy quark light-front Hamiltonian which contains 
explicitly a confinement interaction at long distances and a Coulomb-type 
interaction at short distances. With this light-front QCD Hamiltonian, we 
further demonstrate that the nonperturbative QCD dynamics of the strongly 
interacting heavy hadron bound states can be treated as a weak-coupling 
problem. 
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I. INTRODUCTION 



There are two fundamental features in nonperturbative QCD, the quark confinement and 
the spontaneous breaking of chiral symmetry. These two features are the basis for solving 
hadronic bound states from QCD but none of them has been completely understood. Re- 
cently, Wilson et al. proposed a new approach to determine hadronic bound states from 
nonperturbative QCD on the light-front with a weak-coupling treatment (WCT) The key 
to eliminating necessarily nonperturbative effects is to construct an effective QCD Hamil- 
tonian in which quarks and gluons have nonzero constituent masses rather than the zero 
masses of the current picture. The use of constituent masses cuts off the growth of the 
running coupling constant and makes it conceivable that the running coupling never leaves 
the perturbative domain. The WCT approach potentially reconciles the simplicity of the 
constituent quark model with the complexities of QCD. The penalty for achieving this weak- 
coupling picture is the necessity of formulating the problem in light-front coordinates and 
of dealing with the complexities of renormalization. 

Succinctly, this new approach of achieving a QCD description of hadronic bound states 
can be summarized as follows Using a new renormalization scheme, called similarity 
renormalization group (SRG) scheme 0, one can obtain an effective QCD Hamiltonian H\ 
which is a series of expansion in terms of the QCD coupling constant, where A is a low 
energy scale (~ 1 GeV). Then one may solve from Hx the strongly interacting bound states 
as a weak-coupling problem. The WCT scheme contains the following steps: (i) Compute 
explicitly from SRG the H\ up to the second order and denote it by Hxq as a nonperturbative 
part of H\. The remaining higher order contributions in H\ are considered as a perturbative 
part H\i. (ii) Introduce a constituent picture which allows one to start the hadronic bound 
states with the valence constituent Fock space. The constituent quarks and gluons have 
masses of a few hundreds MeV, and these masses are functions of the scale A that must 
vanish when the effective theory goes back to the high energy region, (iii) Solve hadronic 
bound states with H\q nonperturbatively in the constituent picture and determine the scale 
dependence of the constituent masses and the coupling constant. The couphng constant g 
now becomes an effective one, denoted by g\. If we could show that with a suitable choice 
of A at the hadronic mass scale, the effective coupling constant g\ can be arbitrarily small, 
then WCT could be applied to H\ such that the corrections from H\i can be truly computed 
perturbatively. If everything listed above works well, we may arrive at a weak-coupling QCD 
theory of the strong interaction for hadronic bound states. 

With the idea of SRG and the concept of coupling coherence 0], Perry has shown that 
upon a calculation to the second order, there exists a logarithmic confining potential in the 
resulting light-front QCD effective Hamiltonian [Q. This is a crucial finding to light-front 
nonperturbative QCD. However, the general strategy of solving hadrons through the WCT 
scheme has not been examined. In this paper, I will use SRG to analytically derive from 
the light-front heavy quark effective theory (HQET) a heavy quark QCD Hamiltonian 
which is responsible to heavy hadron bound states. The resulting Hamiltonian explicitly 
contains a confining interaction between heavy quark and antiquark at long distance plus a 
Coulomb-type interaction at short distance. Based on this effective QCD Hamiltonian, I then 
study the strongly interacting heavy hadronic bound states, from which I can qualitatively 
provide a realization of WCT to nonperturbative QCD on the light-front. 
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The reason of choosing heavy hadron systems as a starting example in this investigation 
is that for hght quark systems, both quark confinement and spontaneously chiral symmetry 
breaking play an essential role to the quark dynamics in hadrons. In other words, to provide 
a good QCD description for light quark systems, it is necessary to understand the underlying 
mechanism for quark confinement as well as for chiral symmetry breaking. This will certainly 
complicate the study of nonperturbative QCD. However, for heavy quark systems, chiral 
symmetry is explicitly broken so that confinement is the sole nontrivial feature influencing 
heavy quark dynamics. One may argue that the mass scales for heavy and light quark 
systems are different. The heavy quark energy cannot run down to the usual hadronic scale 
of light quark systems due to heavy quark mass. Meanwhile, confining interactions must 
be energy scale dependent. Apparently, confinement for heavy quark systems could be very 
different from light quark systems. However, despite heavy and light quarks, confinement 
arises only from low energy gluon interactions. In other words, the confinement mechanism 
must be the same for both heavy and light quark systems. We thus choose heavy hadron 
systems without any loss of generality. 

In order to avoid the possible confusion arisen from different mass scales and to cor- 
rectly extract confining interactions in heavy quark dynamics, it is convenient to work with 
heavy quark effective theory (HQET). HQET is a theory of QCD in l/mq expansion [TU 



where tuq is the heavy quark mass. In HQET, the nonperturbative dynamics is determined 
through the interacting gluons and heavy quarks by exchanging a small residual momentum 
of heavy quarks, which is of order Aqcd- As a result, within HQET we can indeed explore 
the nonperturbative QCD dynamics for heavy quark systems in the same scale as that for 
light quark systems. Meanwhile, the extension of the study to light quark systems is straight- 
forward, although undoubtedly the corresponding result must be very complicated due to 
the spin dependence of the nonperturbative interacting Hamiltonian. The spin dependent 
interactions on the light-front are essentially related to the chiral symmetry breaking. These 
spin dependent interactions in HQET are suppressed in the leading order approximation be- 
cause they are l/rriQ corrections and can be treated perturbatively with respect to heavy 
hadron states. This is why for heavy quark systems the chiral symmetry breaking can be 
treated separately from the confinement. 

In fact, the model-based theoretical investigations on heavy quarkonia lasted for one and 
half decades is recently replacing by first-principles exploration on QCD. The lattice QCD 
simulation may give an acceptable description for heavy quarkonium spectroscopy with man- 
ageable control over all the systematic errors 0]. The development of nonrelativistic QCD 
provides a general factorization formula to quarkonium annihilation and production pro- 
cesses so that a rigorous QCD analysis may become possible . Meanwhile, in the past five 
years considerable progress has been made for heavy hadrons containing one heavy quark, 
due mainly to the discovery of the heavy quark symmetry (HQS) ^ and the development 
of HQET |]10[ from QCD. The HQS and HQET have in certain contents put the description 
of heavy hadron physics on a QCD-related and model-independent basis. Moreover, HQET 



has also been extended to describe heavy quarkonia |TT]. Yet, a truly first-principles QCD 
understanding of heavy hadrons is still lacking because so far none is able to give a direct 
computation of heavy hadron bound states from QCD. On the other hand, in the last decade, 
the investigations of the light-front field theory on nonperturbative bound state problems 
have made some progress. Starting with heavy hadrons may provide a possible explicit 
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solution of hadronic bound states in light-front QCD. Hence, in this paper we shall mainly 
concentrate on heavy quark systems, especially for heavy quarkonia, but the generality of 
nonperturbative QCD dynamics obtained in this formalism will also be discussed. 

The paper is organized as follows. In Section 2, the general procedure of constructing 
a renormalized effective Hamiltonian H\ in the light-front SRG is briefly reviewed and the 
possible existence of light-front confining interactions in this formulation is discussed. In 
Section 3, by applying the general procedure to the light-front HQET of QCD, a light-front 
heavy quark confining Hamiltonian is analytically derived. A light-front picture of quark 
confinement is illustrated. In Section 4, the light-front heavy hadronic bound states are 
explored within the WCT scheme. As an example, the light-front heavy quarkonium bound 
state equation is solved by the use of a Gaussian-type wavefunction ansatz, from which the 
WCT to nonperturbative QCD is explicitly explored. Finally, the physical implications in 
the realization of WCT are discussed in Section 5. 



II. EFFECTIVE QCD HAMILTONIAN IN THE SRG SCHEME 

A. Light-front SRG Scheme 

We begin with the general formulation of the similarity renormalization group approach 
to construct a low energy QCD Hamiltonian which was first proposed by Glazek and Wilson 
0. The basic idea of the SRG approach is to develop a sequence of infinitesimal unitary 
transformations that transform an initial bare Hamiltonian to an effective Hamilto- 
nian Hx in a band-diagonal form relative to an arbitrarily chosen energy scale A: 

Hx = SxH^'Sl (2.1) 

Here the band-diagonal form means that the matrix elements of Hx involving energy jumps 
much larger than A will all be zero, while matrix elements involving smaller jumps or two 
nearby energies remain in Hx- The similarity transformation should satisfy the condition 
that for A — >• cx), Hx —>■ H^ and Sx — >■ 1. 

In this paper, we shall follow the formulation of SRG developed on the light-front [|1|. 
The effective Hamiltonian we seek is Hx with A being of order a hadronic mass (~ 1 GeV). 
We begin with a given bare Hamiltonian which can be written by H^ = Hq + H^, where Hq 
is a bare free Hamiltonian and Ei is its eigenvalue. Consider an infinitesimal transformation, 
then Eq.(|2.1j) is reduced to 

(2.2) 

which is subject to the boundary condition lim^^oo-f^A = H^, where Tx is a generator of 
the similarity transformation. 

To force the Hamiltonian Hx becoming a band-diagonal form in energy space, we need 
to specify the action of Tx. This can be done by introducing the scale A with xxij = ^^^\_x 
into a smearing function fxij = f{xxij) such that 

r 1 X < 1/3 

f{x) = I smoothly from 1 to 1/3 < x < 2/3 , (2.3) 
[ X > 2/3 
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and reexpressing Eq.( p.2|) by 

= f\ij[Hi\,Tx\ij + — {In fxij)Hxij, 
Tm, = - M[Hix,Txh - ^(ln/A.,)i^A.,j. (2.4) 

Here we have written Hx = Hq + Hjx because Hq is invariant under transformations. And 
we have also used the notation Aij = {i\A\j), where \i) and \j) are eigenstates of Hq. Since 
f{x) vanishes when x > 2/3, one can see that Hxij does indeed vanish in the far off-diagonal 
region. It also can be seen that Txij is zero in the near-diagonal region. The solutions for 
Hjx and Tx are 

Hix = Hf^ + [Hix'j Ty] , Tx = Hfyj, + [Hiy, Ty] , (2.5) 

where = fxijHf^,, Hf^^^^ = --^^ (^/a.,) Hf^^, and 

Xx'ij = —fxij / dXXyij, (2.6) 
1 / d \ f°° 1 

From eq.( p.5|) , one obtains an iterated solution for Hx, 
Hx = (^Ho + Hf^ + (^ jHfy,Hfyj.] l^ 

+ [[Hfy, Hfy,j,] , Hfyj.] + [Hfy, [Hfy,, Hfy,j,] ] \ + ... 
\ ' V 'Rl " V 'T' 

'R 'R 

= i/f ^ + H^ + i/f + . . . . (2.8) 

Thus, through SRG, we eliminate the interactions between the states well-separated in 
energy and generate the effective Hamiltonian of eq.( |2.8| ). The expansion of eq.( |2.8| ) in terms 
of the interaction coupling constant brings in order by order the full theory corrections to 
this band diagonal low energy Hamiltonian. 

Practically, we do not really compute Hx to all the orders in the expansion of eq.( |2.8| ). 
In the WCT scheme, we consider the first two terms in (|2.8| ) as the nonperturbative part of 
Hx in the determination of hadronic bound states. Thus, in the fist step of WCT, we shall 
only compute Hx up to the second order in coupling constant. We hope that the remaining 
part can be handled perturbatively, as we have mentioned in the introduction and will be 



discussed more later. Up to the second order, eq. (|2.8|) can be further expressed explicitly by 



The front factor in the above equation indicates that Hx only describes long distance in- 
teractions (with respect to the scale A) which is responsible to hadronic bound states. The 
function gxij is given by 
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gxijk 



roo d 



More explicitly, the bare Hamiltonian input in the above formulation can be ob- 
tained from the canonical Lagrangian with a high energy cutoff that removes the usual UV 
divergences. For light-front QCD dynamics, the bare Hamiltonian in our consideration is 
the canonical light-front QCD Hamiltonian that can be either obtained from the canonical 



procedure in the light-front gauge [T^IT^ or generated from the light-front power counting 
rules fl]] . Instead of the cutoff on the field operators which is introduced in ref . 0] , we shall 
use in this paper a vertex cutoff to every vertex in the bare Hamiltonian: 

9{AyP^-\pr-pj\), (2.11) 

where p~ and pj are the initial and final state light-front energies respectively between the 
vertex, A is the UV cutoff parameter, and the total light-front longitudinal momentum 
of the system we are interested in. The theta function is defined as: 6{x) = 1, 1/2 and 
for X > 0, = and < 0, respectively. Eq. (|2.11|) is also called the local cutoff in light-front 
perturbative QCD [|r^. All the A-dependences in the final bare Hamiltonian are removed by 
the counterterms. The use of eq.( |2.11| ) largely simplifies the analysis on the original cutoff 
scheme in ref. |jl|. 

Meanwhile, for a practical SRG calculation, we also have to give an explicit form of the 
smearing function fxij. One of the simplest smearing functions that satisfies the requirements 
of SRG is a theta-function: 

fMj=0{\-XMj). (2.12) 

\p~-p~\ 

On the light-front, it is convenient to redefine xxij = p-_^_p-_^_l2/p+ the light-front. Then 
we can further replace the above smearing function by the following form: 

/a.,=^(;^-|AP,7|), (2.13) 

where AP^y = P^^ — P^ is the light-front free energy difference between the initial and final 
states of the physical processes. The light-front free energies of the initial and final states 
are defined as sums over the light-front free energies of the constituents in the states. 

Note that the choice of a specific smearing function corresponds to picking up a partic- 
ular renormalization scheme to specify the renormalization scale in SRG. This is similar to 
the specification of the renormalization scale in usual RG by picking up a particular regu- 
larization scheme. As we will see later, physics must be independent of the choice of the 
smearing function. With the smearing function given by eq.( p.l3| ), eq.( p.9|) can be simply 
reduced to 



+ 



|AP,^|-AVP+)^(|AP,^|-|APT|) 



APr,\-Xyp^)e{\Ap-\-\APr,\) 



APi, 



API 

+ ■■■). (2.14) 
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B. Application to light-front QCD 



Next, we shall calculate explicitly the low energy effective light-front QCD Hamiltonian 
in the similarity renormalization group scheme. Light-front QCD is the theory of QCD 
defined on the light-front with the light-front gauge. The canonical theory of the light- 
front QCD was formulated in 70's [|15|,|16|. In early 80's, Lepage and Brodsky developed 



the Hamiltonian formulation of light-front QCD and explored extensively the asymptotic 
behavior of hadrons and the perturbative expansion of various exclusive processes . Here 
we shall use the two-component Hamiltonian formulation of the light-front QCD [l^ that 
has been used in illustrating the ideas of the WCT scheme to nonperturbative QCD in the 
previous pubhcation 

The light-front coordinates are defined as follows: = a;° ± x^,a;\ = a;*(z = 1,2), with 
being chosen as the light-front "time" direction along which the states are evolved. Then 
x~ and x± become naturally the longitudinal and transverse coordinates. The inner product 
of any two four-vectors is given by a^a^ = ^{a^b~ + a^h^) — ctj, ■ and the time and space 
derivatives {d^^ = and the 4-dimensional volume element are written as = 2^|^, 

d~ = 2-^, = and d'^x = ^dx^dx~ cPx±. 

In the above coordinate system, the fermion field can be separated hy = ■?/'+ + "4^- 
with = ^'j^'j^iIj, and tp^ and ip^ are the so-called dynamical and constraint components 
of the quark field on the light-front, respectively. Due to the gauge symmetry, only two 
components of the 4-vector gauge field are dynamical field variables. With the choice 
of the light-front gauge A'^ = 0, we can explicitly solve the unphysical components of the 
quark and gauge fields {A~,ip-) from the constraint equations which are obtained from QCD 
Lagrangian Cqcd = -{F^^'Fa^.y + + mq)^): 

= (^)(^ax-Dx + /5m,)^+, ax = 7V, = 7° (2-15) 



A- = 2 



(2.16) 



where Ta is the generator of SU(3) color group: [Tq, T^] = if°-'^'^Tc. Then the QCD Lagrangian 
on the light front can be expressed in terms of only the physical field variables (A^, ■?/'+): 

C = + - (2.17) 

where Ti is the light-front QCD Hamiltonian density [jl3 



n = no + ni, (2.18) 

and 

no = lid'AiWAi) + (2.19) 

= T^ggg + T^ggg + ^-gggq + T^gggg + T^-gggg, (2.20) 

with 



7 



= 9^ {-2 (^) {d-A^) + &. (i^) {a -8^ + 



7^ 



555 



9^ 



+ 2 



9+ 



9^ 



(Aid^Ai;) 



5555 



9 jabc jade ^ 



4* /I* /I-' 



1 



d+ 



+2 - - 



9^ 



(2.21) 
(2.22) 

(2.23) 
(2.24) 

(2.25) 



Here, we denoted A± = T"-Aa_L, Aa± = (A^, A^). The notation a is defined by: = a^, cr^ = 
— (the 2x2 Pauh matrices), and is the two- component form of the hght-front quark 
field: 



■0+ = A"^^ 



V"- = A 







which comes from the use of the hght-front 7-representation: 



(2.26) 



7 



7 



-i 
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, 7^ = 


i 




















, 7^ = 








Thus the bare hght-front QCD Hamihonian is given by 

= J dx~^cfx±(Ho + 'Hpj + counterterms. 



(2.27) 



(2.28) 



where means that the local cutoff, eq. (|2.11| ), has been imposed. Also, counterterms are 
added to eq. ( ^.28] ) in order to remove all the A-dependence. Now, we can see that the 
canonical light-front QCD Hamiltonian is expressed purely in terms of the transverse gauge 
field components A± and the physical two-component quark fields ^. 

Upon the calculation to the second order in the coupling constant, the effective Hamil- 
tonian (|2.14| ) in qq sector (with the initial and final states, \i) = b^{pi, Xi)S{p2, X2)\0) and 
\j) = &^(p3, A3)(i^(p4, A4)|0), respectively, where Pi and Aj denote the respective momentum 
and helicity of a quark on the light-front) becomes 



H 



\ij 



9{\yP^-\APr\)l{J\■.H^■.\^) 



27C 



1 

p+ 



9 

X'^Cf-;^ Ine + mass counterterms 
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XL 



^q]_ ■ P3± - im ^^,a ■ pi± + im 
L — \ — i- o — a - 



q]_ a ■p2±+im^-, ^^la -pi^-im 

I — : 1 — r-; O — O ■ 



\vt\ 



\vt\ 

X-x,Fri\. (2.29) 



In eq.( p.29| ), P+ = V\ ^ V2 = Ps + v\i ^^ij = Pi ~^ P2 ~ Ps ~ Pii '■ : represents a 
normal ordering, where the instantaneous interaction contribution to the quark self-energy 
has been included in the self-energy calculation which is given by the mass counterterms 
and the logarithmic divergence, the color factor Cj = (T"T") = (A^^ — 1)/2N, N = 3 the 
total numbers of colors, and e is an infrared longitudinal momentum cutoff. Since In e is an 
infrared divergence, it cannot be removed by mass counterterms. In gauge symmetry, this 
divergence must be canceled in the physical sector (and this is true as we will see later). The 
last term in ( |2.29[ ) corresponds to one transverse gluon exchange contribution in the second 
order of the low energy Hamiltonian expansion (|2.14|) . The momentum q is carried by the 
exchange gluon: g"*" = pf —pt = pt ~pt: 1± = Pi J- ~P3i. = P4± ~P2L, XXi denotes a hehcity 
eigenstate. The factor Frij arises from the similarity renormalization transformation: 



eiA^p^ - K -P3 - q-m^VP^ - \P,-P2 - q-\) 



0i\Pi -Ps -Q I - pt)^(IPi -P3 -Q \-\Pi -P2 -Q I) 

Pi -P3 - 

, 0{\P4 -P2 - - jrWilPi -P2 - Q'l - \Pi -P3 - Q'l) 



Pi -P2 -Q 

+ [pi ' — ^ P3 , P2 < — > Pi , q — > -q) 



(2.30) 



Pi 



and q 
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All light-front energies in eq.( p.30| ) are on mass-shell: p^ 

If we continue to evaluate all the terms in the expansion of eq.(|2.14|), the resulting 



Hamiltonian is the exact QCD Hamiltonian at the energy scale A. In practice, we only 
consider the leading and the next-to- leading terms, i.e., eq.( p.29| ), as a starting effective 
Hamiltonian. The basic idea to realize a weak-coupling treatment of QCD for hadrons is 
whether we can solve hadron states from this effective Hamiltonian (|2.29|) with an arbitrary 



small coupling constant g such that the higher order corrections in (|2.14| ) can be handled 
perturbatively. From the success of the constituent quark model, we understand that a 
necessity for such a realization might depend on the existence of a confining interaction in 

(Pi). 



C. Quark confining interaction on the liglit-front 

Naively, we know that any weak-coupling Hamiltonian derived from QCD will have 
only Coulomb-like interactions, and confinement can only be exhibited in a strong-coupling 
theory. However, Perry has recently found that upon to the second order calculation of 
the low energy Hamiltonian, a logarithmic confining potential has already occurred 0. 
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Explicitly, the two-body quark- ant iquark interaction from the first term of eq. (|2.29|) is the 
instantaneous gluon exchange interaction Tigggg which has the form in the momentum space: 

(2.31) 



The confinement must be associated with the interaction where g+ 0. This is because 
only these particles with zero longitudinal momentum can occupy the light-front vacuum 
state |l[]. Thus, we need to analyze the feature of the effective Hamiltonian when — 0. 
For 0, the dominant second order contribution from one transverse gluon exchange 

interaction in eq.( p.29| ) is given by 

1 ll , 1 



. +^2 2 — ] 



e{q--y/P+). (2.32) 



In the usual perturbative calculation, no such a ^-function (coming from the smearing func- 
tion) is attached in the above equation. Thus these two dominant contributions from the 
instantaneous and one-gluon exchange interactions are exactly cancelled when g"*" — >■ 0. Only 
Coulomb-type interaction (the terms ~ 0(^)) remains. However, in the light-front sim- 
ilarity renormalization group scheme, the corresponding one-gluon exchange contribution 
of eq.( p.32 ) only contains these gluons with energy being greater than the energy cutoff 



\^ I p+ _ As a result, the instantaneous gluon exchange term l/(g"'")^ remains uncancelled 
if the gluon energy, q]_/q^., is less than A^/P"*". The remaining uncancelled instantaneous 
interaction contains an infrared divergence and a finite part contribution (for a detailed 
derivation, see the next section). The divergence part is cancelled precisely for physical 
states by the same divergence in the quark self-energy correction [see ( 2.29|) ]. The remain- 
ing finite part corresponds to a logarithmic confining potential: 

6Aln(P+|x-|) + CAln(AV±n- 

The above result is indeed first obtained by Perry with the use of the concept of coupling 
coherence and a slightly different renormalization scheme (also see a oversimple derivation 
given by Wilson [0). Here the derivation is purely based on the light-front similarity 
renormalization scheme. 

One may argue that the existence of such a confining potential in H\ up to the second 
order may only be an artificial effect designed in the renormalization scheme we used here. 
If we included the interaction with the exchange gluon energy below the cutoff, then the 
instantaneous interaction would be completely cancelled, and no such confining potential 
should exist, as expected in the usual perturbation computation. Wilson has pointed out 
that the set up of the new renormalization scheme is motivated by the idea that the gluon 
mass must be nonzero in the low energy domain (a constituent picture). The massive gluon 
can be originated from the nonlinear interactions in non-abelian gauge theory. Therefore, 
the non-cancellation of the instantaneous interaction in the low energy domain is indeed a 
consequence of the existence of the constituent massive gluons due to the non-abelian gauge 
interactions. It is independent of any particular renormalization scheme. The use of the low 
energy cutoff A (namely, taking the smearing function as a theta function) just gives us a 
simple realization of this confining picture that the exchange energy of massive dynamical 
gluons cannot run down to the zero value in nonperturbative QCD. 
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Yet, although there exists a confining interaction in Hx even up to the second order, 
this effective QCD Hamiltonian is aheady very comphcated, due to the spin dependent 
part in the interactions. As we know the spin dependent interactions on the hght-front are 
essentially related to the chiral symmetry breaking. In order to examine the above picture of 
confining mechanism on the light-front and to develop explicitly a weak-coupling treatment 
approach of nonperturbative QCD to hadronic bound states, in the next section we shall 
utilize the above formulation to heavy quark systems. We find that the low-energy QCD 
Hamiltonian Hx for heavy quarks can be largely simplified and an analytic form consisting 
of the confining and Coulomb potentials emerges. 



III. HEAVY QUARK CONFINING HAMILTONIAN 

In the past few years, QCD has been made a numerous progresses in understanding the 
heavy hadron structure, due mainly to the discovery of heavy quark symmetry by Isgur 
and Wise 0, and the development of HQET by Georgi et al. [l^. However, to completely 



understand the nonperturbative QCD dynamics of heavy hadrons, one should solve heavy 
hadron bound states directly from QCD, which is still an open question. Next we will use 
SRG to the light-front HQET to derive a heavy quark confining Hamiltonian, from which 
we may solve from QCD the heavy hadron bound states directly. 

A. Light-front HQET 

The light-front heavy quark effective Lagrangian derived from QCD lagrangian C = 
Q{ip — rnQ)Q as a l/mg expansion is given in refs. [^||: 

9 oo -| 

C = —QUzv ■ D)Q,+ - ( ^)"qU{(^« ■ D){-zD+r-\za ■ (3.1) 

n=l ''^OV 



where Qy+ is the light-front dynamical component of the heavy quark field after the phase 
redefinition: Q{x) = e~^"^'^'"'^ {Qv+{x) + Qv~{x)), the four velocity of the heavy hadrons, 
= MffV^ with = 1 and the heavy hadron mass, a-D = Q:±-D_i — ^(a±-v±+P)D^, 
and D^^ the usual covariant derivative. The corresponding light-front heavy quark bare 
Hamiltonian density is given by 

n = -^QU{v-d+ - 2v^ ■ d^)Q,+ - ^QU{v ■ A)Q,+ 

oo 1 

+ E ( ■ D){-tD^r~\ta ■ D)}Q,+ . (3.2) 

n=l "^Q"^ 

The heavy antiquark Hamiltonian has the same form except for replacing v by —v. 

In the large mg limit, only the leading (spin and mass independent) Hamiltonian is 
remained. The l/mq terms {n > 1) in ( p.2|) can be regarded as perturbative corrections to 
the leading order operators and states. To determine confining interactions in heavy quark 
systems, the leading heavy quark Hamiltonian plays an essential role. With the light-front 
gauge A'^ = 0, the leading-order bare QCD Hamiltonian density is 
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1 



9^ 



+ (^) (3-3) 

where '?/^+ is either the heavy antiquark field or the hght-front quark field operator in the 
present consideration. Note that besides the leading term in eq. (|3.2| ), the above bare Hamil- 
tonian has also already included the relevant terms from the gauge field part, — iTr(F^j,F'^'^), 
of the QCD Lagrangian. These terms come from the elimination of the unphysical gauge 
degrees of freedom, the longitudinal component A~ |]13[. Eq. ( |3.3| ) has obviously the spin 
and fiavour heavy quark symmetry, or simply the heavy quark symmetry. 

The above leading Hamiltonian (or Lagrangian) is the basis of the QCD-based description 
for heavy hadrons containing a single heavy quark, such as B and D mesons. As recently 
pointed out by Mannel et al. [0,0 the purely heavy quark leading Lagrangian may be not 
appropriate to describe heavy quarkonia. This is because the anomalous dimension of QCD 
radiative correction to QQ currents contains an infrared singularity in the limit of two heavy 
constituents having equal velocity. Such an infrared singularity is a long distance effect and 
should be absorbed into quarkonium states. To solve this problem, they show that one 
should incorporate the effective Hamiltonian with the first order kinetic energy term into 



the leading Hamiltonian [11|. The light-front kinetic energy can be obtained from eq. (|3.2| 



n,,n = — ^—QlJdl - ^-^^d^ + ^9+4q.+. (3.4) 

Indeed, as we will see in Section IV. C, for heavy quarkonia the heavy quark carries a relative 
momentum that is proportional to heavy quark mass. Then the first order kinetic energy 
term is of the same order as the leading Hamiltonian. Other terms in the l/mq are really 
suppressed by the power of I/itlq. This is why one have to incorporate the leading Hamil- 
tonian with the first order kinetic energy term. As a consequence, in the heavy mass limit, 
quarkonia have spin symmetry but no fiavour symmetry. 



B. Confining Hamiltonian for heavy quarkonia 

Within light-front HQET, we now follow the procedure described in the previous section 
to find an effective QCD Hamiltonian for QQ systems. The bare Hamiltonian for QQ systems 
contains ( |3.3| ) and ( p.4| ) for both heavy quark and antiquark plus the full QCD Hamiltonian 
for gluons and light quarks [|l^]. Since the kinetic energy ( p.4|) is of order KqcdI'^q-, which 



is at most of the same order of the Coulomb interaction, we can treat the kinetic energy in 
the same way as the instantaneous QQ interaction [the last term in eq. (|3.3| )]. Thus, the free 
Hamiltonian Hq used in SRC is given only by the first term in eq. (|3.3| ) plus the free gluon 
Hamiltonian. 

Keeping the above consideration in mind, it is easy to compute the effective Hamiltonian 
eq.( |2.14| ) for QQ systems. Following the WCT ideas, we shall calculate Hx for QQ systems 
up to the second order in the initial and final states defined by |i) = h\{k\^ \\)d\{k2-, -^2)|0) 
and |j) = h\{k-^,\-^d\{k^,\/^\^) , respectively, where fcj is the residual momentum of heavy 
quarks, = mqv^ + fcf, and Aj its helicity. The operator h\{ki^\i) [(i^(A;j, Aj)] creates a 
heavy quark [antiquark] with velocity f , 
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{K{k, A), blik', A')} = A), dl,,{k', A')} = 2{2n)%,,5xx'5%k - F), 

and _ k') = 6{k+ - k'^)5'^{k^ - k'^). The result is 

where 



^QQfreeij + ^QQij^ 



(3.5) 



(3.6) 



^QQfreeif- 



[2{27ryY6'{k, - k,)6'{k, - k^)6x,x,6x,x. 



X 



A r 



2Ki + A^{2y'-2y + l] 



Vq-Qijiy - y' , 



2{2^f5''{h + k2-h 



A' - 2^C, ^ 



In e 



^4)'^AiA3^A2A4 



X 



-Ag^{T''){T''' 



+ 



{K^-K'^r + iy-y'^A 



eiAiy-y\ - k'^, A) ~ X') 



(3.7) 



(3.8) 



Here we have introduced the longitudinal residual momentum fractions and the relative 
transverse residual momenta, 



y 
y' 



k+/K+ 



fi± = ki^ - yK^, 
«j_ = k3± - y'Kj 



(3.9) 



where K^^ is defined as the residual center mass momentum of the heavy quarkonia: K'^ = 
Av^, and A = Mh — faq — ttiq is a residual heavy hadron mass. It follows that = 
ki + k2 = k^ + k^, K± = ki± + k2± = k^^ + ki±. Since Q < pt = ttt^qV^ + kf < Mhv^, 
in the heavy quark mass limit, we have Mh — > 2mQ so that —rnqv^ < k^ , kf < rnqv^ . 
Hence, the range of y and y' are given by — oo < y,y' < oo. We have also defined in eq, 



A{y - y' , K± - k'^, A) = 



K 



\y-y'\ 



+ \y-yW 



(3.10) 



Eq. (|3.6| ) is the nonperturbative part of the effective Hamiltonian for heavy quarkonia 
in the WCT scheme, in which we have already let UV cutoff parameter A — *■ cxd and the 
associated divergence has been put in the mass correction. The kinetic energy (|3.4| ) is now 
included in the above effective Hamiltonian [the l/mq term in eq.( |3.7] )]. Note that there is 
an infrared divergent term in eq.( p.7|) which comes from the quark self-energy correction in 
SR.G, 



2S 



2tx' 



25ml, 



(3.11) 



where e is an infrared cutoff of the momentum fraction / K'^ , and the longitudinal 
momentum carried by gluon in the quark self-energy loop. The usual mass correction dniq = 



^CfA In^, has been renormalized away in eq. (|3.7| ), where the color factor Cf = (T°'T"-) 



(N —1)/2N, N = 3 the total numbers of colors. In the WCT scheme, by removing away this 
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mass correction, we should assign the corresponding constituent quark mass in Hxo being 
A-dependent. But, the heavy quark mass is larger than the low energy scale. Its dependence 
on A should be very weak and could be neglected. Meanwhile, the QQ interaction ( p.8| ) 
contains two contributions: the instantaneous interaction plus the second order contribution 
in eq.( p.l4| ) [i.e. the terms proportional to the theta function in eq. (|3.8| )]. We shall show next 
that the above Vqq is indeed a combination of a confining interaction plus a Coulomb-type 
interaction. 



C. A simple picture of quark confinement on the light-front 

In the conventional picture, QCD has a complex vacuum that contains infinite quark 
pairs and gluons necessary for confinement and chiral symmetry breaking. On the light- 
front, the longitudinal momentum of physical particles is always positive, = > 0. 
Consequently, only these constituents with zero longitudinal momentum (called zero modes 
pop can occupy the light-front vacuum. The zero modes carry an extremely high light-front 



energy which has been integrated out in SRG. As a result, some equivalent effective interac- 
tions are generated in H\ so that the light-front QCD vacuum becomes trivial. The nature of 
nontrivial QCD vacuum structure, the confinement as well as the chiral symmetry breaking, 
is then made manifestly in Hx in terms of these new effective interactions. We will see that 
HxQ explicitly contains a confining interaction at long distances. The interactions associated 
with the chiral symmetry breaking may be manifested in the fourth order computation of 



Hx for light quark systems |[T^, but these interactions are not important in the study of 
heavy hadrons here. 

The confining interaction can be easily obtained by applying the Fourier transformation 
to the first term in ( |3.8| ). It is convenient to perform the calculation in the frame Kj_ = 0, 
in which 



2 

9l f^^'' ^^+^iq+x- llm 2Jl(|Xx|gXm) ^3^2) 



Stt^ Jo \x±\q±m 
where = kf — = K^{y — y'),qi_ = ki± — ks_L = k,± — for Kj_ = 0, q± 



7-2 



'W+'i^ ~ 'W+^'l^'^' ^^"^ is a Bessel function. An analytic solution to the integral (|3.12|) 
may be difficult to carry out. However, the nature of confining interactions is a large distance 
QCD behavior. We may consider the integral for large x~ and x±. In this case, if q'^x~ 
and/or |x_L|g_Lm are large, the integration vanishes, yet Ji(x) = | + jg + ■ ■ ■ for small x. The 
dominant contribution of the integral (|3.12|) for large x~ and x±_ comes from the small g"*" 
such that q'^x~ and/or |x_L|g_Lm, must remain small, which leads to e^'^^^ 2Ji(|a:x|g±m) ^ -^^ 

This corresponds to < ^ and/or g+ < 



If g+ < ^ < , j2;,2 , eq. (|3.12| ) is reduced to 




ln(is:+|x-|) + lne), (3.13) 
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where a term ~ ^ is neglected since x~ is large, and e is an infrared cutoff of the momen- 
tum fraction / . It is the same as in eq. (|3.11|) so that the above infrared logarithmic 
divergence (~ Ine) exactly cancels the divergence in eq.( p.7|) for color single states. What 
remains is a logarithmic confining interaction (except for a color factor): 



2 \2 



V.onf{x-,x^) ~ 7^^HK'-\x~\). (3.14) 



Similarly, when g"*" < ^J^^y^2 < have 

g2 . 1 . \2 —2 „2\2 



2,4"^ "'^^^(1^'^ - 7^'") - ^(MAVxlVln.). (3.15) 

where the term ~ has also been ignored because of the large x\. Again, the infrared di- 
vergence (~ Ine) is cancelled in Hx for physical states, and we obtain the following confining 
interaction: 

V.onf.{x-,x^) ~ ^fl^MA^I^J^). (3.16) 

Hence, the effective Hamiltonian Hxo exhibits a logarithmic confining interaction between a 
heavy quark and a heavy antiquark in all the directions of x~ and x± space. 

The Coulomb interaction corresponds the second term in ( |3.8|) , its Fourier transformation 
(except for the color factor) is 



1 




where . f i + (^)%-)^ which is defined as a "radial" variable in the light-front space 
|l[]. Eq.( p.l7| ) shows that the Coulomb interaction on the light-front has the form 

W(x-,xx)~-|^^-. (3.18) 
47r A+ ri 

Thus, we have explicitly shown that H\q contains a Coulomb interaction at short distances 
and a confining interaction at long distances. 

Moreover, a clear light-front picture of quark confinement emerges here. To be specific, 
we define quark confinement as follows: i) There is a confining interaction between quarks 
such that quarks cannot be well-separated; ii) No color non-singlet bound states exist in 
nature, only color singlet states with finite masses can be produced and observed; and iii) 
The conclusions of i-ii) are only true for QCD but not for QED. 

We have shown explicitly the existence of a confining interaction in H\q. One can also 
easily see from Hxq the non-existence of color non-singlet bound states. This is essentially 
related to the infrared divergences in Hxo- From eqs. (|3.13|) and ( 3.15| ), we find that the 



uncancelled instantaneous interaction contains a logarithmic infrared divergence. Except 
for the color factor, this infrared divergence has the same form as the divergence in eq.(|3T 
Thus, we immediately obtain the following conclusions. 
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(a) . For a single (constituent) quark state, the interaction part of Hxq does not contribute 
to its energy. The remaining infrared divergence from quark self-energy correction implies 
that the dynamical quark mass for a single quark state is infinite (infrared divergent) and 
cannot be renormalized away in the spirit of gauge invariance. Equivalently speaking, single 
quark states carry an infinitely large mass and therefore they cannot be produced. 

(b) . For color non-singlet composite states, the color factor iT"-)ai3(T"')s'y in the QQ inter- 
action is different from the color factor Cf = Tr (T"T"). Therefore, the infrared divergence 
in the self-energy correction also cannot be cancelled by the corresponding divergence from 
the uncancelled instantaneous interaction. As a result, color non-singlet composite states 
are infinitely heavy that they cannot be produced as well. 

(c) . For color singlet QQ states, the color factor (T")(T") — > Cj. Thus, the infrared 
divergences are completely cancelled and the resulting effective Hamiltonian is finite. In 
other words, only color singlet composite are physically observable. 

Finally, we argue that the above mechanism of quark confinement is indeed only true 
for QCD. As we have seen the light-front confinement interaction is just an effect of the 
non-cancellation between instantaneous interaction and one transverse gluon interaction 
generated in SRG. Such a non-cancellation arises in SRG because we introduce the energy 
scale A through the smearing function. Introducing the energy scale A in SRG forces the 
transverse gluon energy involved in the QQ effective interaction never be less than a certain 
value (the energy scale A) . This implies that the gluon may become massive at the hadronic 
mass scale. Of course, such a gluon mass must be a dynamical mass generated from the 
highly nonlinear gluon interactions. In other words, the above confining picture is indeed 
a dynamical consequence of non-Abelian gauge theory. This confinement mechanism is not 
valid in QED. In QED, since photon mass is always zero, the photon energy covers the entire 
range from zero to infinity. Thus, in QED, we can always choose the energy scale A being 
zero. With A = 0, the infrared divergences do not occur in the electron self-energy correction. 
As a result, the renormalized single electron mass is finite, in contrast to the divergent mass 
of single quark states. For the same reason, with A = 0, the instantaneous interaction in 
the effective QED Hamiltonian is also exactly cancelled by the same interaction from one 
transverse photon exchange so that only one photon exchange Coulomb interaction remains. 
Thus, applying SRG to QED and let A = in the end of procedure, we obtain a conventional 
QED Hamiltonian which only contains the Coulomb interaction. 



D. Generality of the confining Hamiltonian in SRG 

As we have pointed out in the previous section, the existence of quark confining inter- 
action in this formalism of nonperturbative QCD should be independent of the choice of a 
particular renormalization scale. In this subsection, we shall examine if such a light-front 
confining picture could be an artificial effect designed in our SRG renormalization scheme. 



especially if it is obtained from our specific choice of the smearing function (|2.13| ). To an- 
swer this question, we shall rederive H\ without specifying the detailed form of the smearing 
function fxij. 

For heavy quarkonia in the HQET, AKr^ = AKr^ = - ^+(^_^,) ((/^± - + (y-y'^A^) ■ 
This immediately leads to fxik = fxjk because xxik = xxjk = k7+k-^^k-^+>:'/k+ - 
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eq. (|2.10|) is reduced to 

gxijk = 1 - fxik = 1 - fxjk- (3.19) 
Following the same calculation of Sec. III. B, we can find that 

Hxoij = fxij{HQQj^^^^. + Vq^.j), (3.20) 

where 



■ A r 



2Ki + A\2y' - 2y + 1)\ -A'-2f-^Cfj^fJnel (3.21) 



^QQijiy - y'^ - i^'a) = 2(27r) 5 {ki + k2 - k^ - ki)6x,XsSx2Xi 

-Ag\T^)iT^) \ 1 A' _ 1 

{K^y I iv - y? ^ {.^ - n'^y + {y- yn'^ ^ ^ ^ 

In Eq.( |3.21| ), fg depends on the detailed form of fxik- 

If we took the smearing function fxik by the theta function of Eg. ( p. 13 ), then Eg. ( p. 22 ) 
will immediately be reduced to Eq.( |3.8| ). Now we shall consider the general case. Note that 

+K~ = 2Av~ « X'^/K'^ since A is of the same order Aqcd — 0.1 ~ 0.4 GeV (for heavy 
quarkonia, A is indeed about 0.05 ~ 0.20 GeV) but A is of the order 1 GeV. Meanwhile, for 
confining interaction, we shall only concentrate in the region = {y — y')K^ — >■ 0. Thus, 

. ^. - - - ^3 23) 

where q± and g+ are the transverse and longitudinal momentum carried by the exchanged 
gluons (see after Eg. p.l2| )). According to the properties of fxij, 

i ql < A2g+/2K+ xxik <\^fx^k = l 
A2g+/2fr+ < gi < 2A^g+/ir+ \ < ^xik < I ^ fxik smoothly goes from 1 to . 



[ ql > 2X'q+/K+ ^ x^k > | - /. 



Xik 



(3.24) 

Then the dominant contribution of the Fourier transformation of Vq-q- for g+ ^ is given 
by 



(27r)2 \ K + 2 ^y_y. 

9x f j„+An+x-(l±i 2Ji(|x^|g^i 



dq+e^^-^- ^"^'7'""^^ + ■ ■ ■ , (3.25) 
271'^ J q+^ \x_L\q_Li 

where g^^^ = X'^q'^ /2K~^ , the dots denote the contribution from the region of | < xxik < §■ 
Introducing a smooth function in this region is motivated for an easy control of possible 
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singularities in the further numerical computations [|T|. However, we found that the use of 
( p.l3| ) for fxik does not cause any new singularity for heavy quarkonium systems. It is more 
interesting to see that the dominant contribution of Vqq^j (the first term in eq. (|3.25| )) has the 
same form of eq.( p.l2| ) and results in the same confining interactions of eqs. (|3.14|) and (|3.16|) 
(except for a factor 1/V2 = fg) that are obtained by using eq.( |2.13| ). This shows that the 
existence of confining interactions in SRG is a general feature in our formulation. The choice 
of eq.( |2.13| ) is just a simple way to introduce a convenient low-energy renormalization scale 
to nonperturbative QCD that makes the confining interaction explicit for heavy quarkonium 
systems. Indeed, Perry and Wilson have used different ways to show the existence of such a 
confining interactions PJT8|] . The crucial point is that the above confining picture is indeed 
model independent since the parameter X is a renormalization scale. Under the consideration 
of SRG (see Section V), physical observables are independent of X. In other words, the 
confining interactions obtained here are derived from QCD without introducing new free 
parameters in SRG scheme. 



E. Extension to heavy-light quark systems 

The heavy-light quark system (heavy hadrons containing one heavy quark) is one of the 
most interesting topics in the current study of heavy hadron physics. We now apply SRG 
to such systems. 

The bare cutoff Hamiltonian we begin with for heavy-light quark systems is the combi- 
nation of the heavy quark effective Hamiltonian (|3.2| ) and the full Hamiltonian for the light 
quarks and gluons |13|. We may also introduce the residual center mass momentum for 
heavy-hght systems, = Af + = pf + kt = pt + k2 , K± = Av± = pi± + ki± = p2i_ + /c2±, 
where A = — mq, pi and p2 are the light antiquark momenta and ki and ^2 the residual 
momenta of the heavy quarks in the initial and final Qq states respectively. 

Following the general procedure, it is easy to find the nonperturbative part of the effective 
Hamiltonian for heavy-light quark systems. 



K+ 



Qqij ( ; 



where 



(3.26) 



HQqfreeij= [2{27ryY S"" {ki - k2)S%Pi - P2)Sx,x,Sx,X, 

^Qqijiv - y'^ '^^ ~ '^J-) = 2{2nf6^{ki +pi-k2- P2)5\^\z^\2\a 
{ 2 



-2^2(^T")(T"^ 



iy - y') 



Av_|_ ^ J_ ) ^_L '^J- ' ^_L ^_L ' ^_L )_L 



{y - y') 



e{B - x^)e{B - A) 



1{K^ - K'^y - {y - y')C-& 



y{y - y') 



+ 



y'{y - y') 

9{A- X^)9{A- B) 



[K^-n'^y + {y-y'yA^\ 



y y 



(3.27) 



(3.28) 
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with B 



and the function A has the same form as in quarkonium 



y-y 

case. Here we have also introduced y = Pi /K^, k± = pi± — yK±, but the range of y is now 

given by < y = < oo. 

The heavy-hght quark effective Hamiltonian is mg-independent. This is because in 
heavy-hght quark systems the heavy quark kinetic energy can be treated as a perturbative 
correction to Hxq. Obviously the above Hxo has the heavy quark spin and flavour symmetry. 
Compared to the Vqq, Vq-q interactions are much more complicated. But it is not difficult to 
check that the above Vq-q contains a confining interaction when is small. The confining 
mechanism is the same for QQ and Qq systems, as well as for qq systems. 

In conclusion, we have obtained in this section the nonperturbative part of a confining 
QCD Hamiltonian for heavy-heavy and heavy-light quark systems in the WCT scheme. We 
have also shown that the confining mechanism is independent of a specific choice of the 
smearing functions. It is indeed also independent of the HQET, namely the same confining 
picture occurs when we include all the l/mq corrections (see the discussion in the previous 
section). We are now ready to study hadron states on the light-front and to explore how 
the WCT scheme works in the present formulation 



IV. THE WCT TO HEAVY HADRON BOUND STATES 

As we mentioned in the Introduction, the ideas of WCT to nonperturbative QCD is to 
begin with the effective QCD Hamiltonian Hx = Hxo + Hxi- Then using the constituent pic- 
ture to solve nonperturbatively the hadronic bound state equations governed by Hxo and to 
determine the running coupling constant gx- If one could properly choose the nonperturba- 
tive Hxo such that gx is arbitrarily small, then the corrections from Hxi could be computed 
perturbatively, and we would say that a WCT to nonperturbative QCD is realized. In this 
section, we shall study such a WCT to heavy hadron bound states. 



A. General structure of light-front bound state equations 

In general, a hadronic bound state on the light-front can be expanded in the Fock space 
composed of states with definite number of particles [|T^,^. Formally, it can be expressed 
as follows 

|vI/(P+,P^,A,)) = ^ I {l[[d'p])2{27if6%P-J2p,)\n,p,X,)<^^{x„K^„X,), (4.1) 

n,Xi i i 

where P~^, P± are its total longitudinal and transverse momenta respectively and its 
total helicity, |n,p. A,) is a Fock state consisting of n constituents, each of which carries 
momentum pi and helicity Aj (J2i K = ^s)] ^{^i, hi±i, A,) the corresponding amplitude which 
only depends on the helicities A,, the longitudinal momentum fractions Xi, and the relative 

transverse momenta Xj = K,i± = pi± — XiP^_. 

The eigenstate equation that the wave functions obey on the light-front is given by 

Hlf\P^, Px, a.) = |P+, P^, A.), (4.2) 
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where Hlp = P the hght-front Hamiltonian. Exphcitly, for a meson wave function, the 
corresponding hght-front bound state equation is: 



- M^) 



{qq\Hint\qq) {qq\H^nt\qqg) 
{qqg\Hint\qq) 



99 



999 



(4.3) 



where Mq 



— the so-called invariant mass, Hint the interaction part of Hip. 



Obviously, solving eg. ( [4 .31) from QCD with the entire Fock space is impossible. A basic 
motivation of introducing the WCT scheme is to simplify the complexities in solving the 
above equation. In the present framework, Hlf = -f^A, where H\ has already decoupled 
from high energy states. Furthermore, the reseparation Hx = Hxo + Hxj is another crucial 
step in WCT, where only Hxo is assumed to have the nonperturbative contribution to bound 
states through eq. (|4.3|) , and Hxi is supposed to be a perturbative term which should not be 
considered when we try to solve eq. (|4.3| ) nonperturbatively. 

The next important step in the WCT scheme is the use of a constituent picture. The 
success of the constituent quark model suggests that we may only consider the valence quark 
Fock space in determining the hadronic bound states from Hxo- In this picture, quarks and 
gluons must have constituent masses. This constituent picture can naturally be realized on 
the light-front However, an essential difference from the phenomenological constituent 
quark model description is that the constituent masses introduced here are A dependent. 
The scale dependence of constituent masses (as well as the effective coupling constant) is 
determined by solving the bound states equation and fitting the physical quantities with 
experimental data. But for heavy quark mess, this A-dependence can be ignored. Once 
the constituent picture is introduced, we can truncate the general expression of the light- 
front bound states to only including the valence quark Fock space. The higher Fock space 
contributions can be recovered as a perturbative correction through Hxi- Thus, eq.( [4.1| ) for 
heavy quarkonia can be approximately written as: 



X(l^oo(y^ f^±Pviki, Xi)dl^{k2, A2)|0), 



A1A2 ' 



(4.4) 



where the wavefunction (pQ-Q^y, k_l) may be mass dependent due to the kinetic energy in 
Hxo [see ( ^.7| )] but it is spin independent in heavy mass limit. Also note that the heavy 
quarkonium states in heavy mass limit are labelled by the residual center mass momentum 



K^. We may normalize eq.( |4.4|) as follows: 

(*(ir'+, i^i, X:)\^{K+, K^, A.)) = 2{27rfK~'6'{K - K')6x'^x. 
which leads to 

dyd'^K± 2 1 



(4.5) 



(4.6) 



With the above analysis on the quarkonium states, it is easy to derive the corresponding 
bound state equation. Let H^p = Hxo of eq.(^]6|), eq. ([4.3|) is reduced to 
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27r2 



+ 



,e{A-\')\ct>^M,K'^). (4.7) 



This is the hght-front bound state equation for heavy quarkonia in the WCT scheme. 
For the heavy mesons containing one heavy quark, similar consideration leads to 



|2A^ - A 



yk- 



4 + m2(A) 



yk 



>^Qq{y,k±, Xl,X2) 



U2 „./ 



27T 



dy'd^fi'j 



■Vqqiy - y', «± - K.'±)^Qq{y', Xi, x 



2), 



(4.8) 



where Vqg is given by eq. 



Note that the light antiquark here is a brown muck, a 



current light antiquark surround by infinite gluons and qq pairs that results in a constituent 
quark mass which is a function of A. 



B. A general analysis of light-front wavefunctions 



A numerical computation to the bound state equations, eqs. (|4^ ) and (|4.8|) , is actually 
not too difficult. However, to have a deeper insight about the internal structure of light-front 
bound states and to determine the scale dependence of the effective coupling constant in 
Hx, it is better to have an analytic analysis. For this propose, we would like to present a 
general analysis of light-front hadronic wavefunctions and then use variational approach to 
solve the bound state equations. 

The heavy hadronic wavefunctions in the heavy mass limit are rather simple. First of all, 
the heavy quark kinematics have already added some constraints on the general form of the 
light-front wavefunction 0(x, k±). When we introduce the residual longitudinal momentum 
fraction y for heavy quarks, the longitudinal momentum fraction dependence in (j) is quite 
different for the heavy-heavy, heavy-light and light-light mesons. 

For the light-light mesons, such as pions, rhos, kaons etc., the wavefunction (f)qq{x, k_i_) 
must vanish at the endpoint x = or 1. This can be seen from the kinetic energy term in 

eq. ([4.3|) , where Mq = — the valence Fock space. To ensure that the bound 

state equation is well defined in the entire range of momentum space, \(pgg{x, must fall 
down to zero in the longitudinal direction not slower than 1/x and 1/(1 — 2;) when a; — >■ 



and 1, respectively. In other words, at least 4>qq{x, k±) ~ \lx{l — x) . For heavy-light quark 
mesons, namely the B and D mesons, the wavefunction (pQqiy, k_l) is required to vanish at 
y = 0, where y is the residual longitudinal momentum fraction carried by the light quark. 
This is because the kinetic energy in eq.( [4.8|) only contains a singularity at y = 0. On the 
other hand, since < y < oo, (t>Qg{y, n±) should also vanish when y ^ oo. Hence, a possible 
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simple solution is (f)Qq{y, k,±) ~ ^y^""^ \/y^~"^ ■ Other form to suppress the singularity 
at y = 0, like e~"/^ instead of ^y, is also possible. For heavy quarkonia, —oo<y<oo, 
the normalization forces (^QQ^y, to vanish as y ^ ±00. Thus, a simple solution may be 

^QQ{y,i^±) ~ e-'^y'. 

On the other hand, the transverse momentum dependence in these light-front wavefunc- 
tions should be more or less similar. They all vanish at k,± ±00. A simple form of the 
K± dependence for these wavefunctions is a Gaussian function: e^'^^^'^^ . 

The above analysis of light-front wavefunctions is only based on the kinetic energy prop- 
erties of the constituents. Currently, many investigations on the hadronic structures use 
phenomenological light-front wavefunctions. One of such wavefunctions that has been widely 
used in the study of heavy hadron structure is the BSW wavefunction , 



4'BSw{x, K,±) = U\Jx{l - x) exp 



2cj2 



exp 



Ml, 



[X 



Xq) 



(4.9) 



2Mi 



where A/" is a normalization constant, u a parameter of order Aqcd, a^o = (| 
and Mh, ttli, and 1712 are the hadron, quark, and antiquark masses respectively. In the 
heavy mass limit, the BSW wavefunction can be produced from our analysis based on the 
light-front bound state equations. 

Explicitly, for heavy-light quark systems, such as the B and D mesons, one can easily 
find that in the heavy mass limit, mi = tuq ~ M^, fng << tuq so that Xq = . Meanwhile, 

we also have Mhx = Ay. Furthermore, the factor \Jx{l — x) can be rewritten by in 
according to the corresponding bound state equation discussed above. Thus, the BSW 
wavefunction is reduced to 



0Qg(?/, k_l) = ^^^/y exp 



2cj2 



exp 



2u' 



-,y 



(4.10) 



This agrees with our qualitative analysis given above. Indeed, using such a wavefunction we 
have already computed the universal Isgur-Wise function in _B — D,D* decays 0: C,{v-v') = 
and from which we obtained the slope of ^(f -f ') at the zero-recoil point, = —^'(1) = 1, 
in excellent agreement with the recent CLCO result |^ of = i.Ql ± 0.15 ± 0.09. 

For heavy quarkonia, such as the bb and cc states, mi = m2 = mg which leads to 
Xq = 1/2 in eq.( 



). Also note that Mh{x — 1/2) = Ay, and the factor \Jx{l — x) must be 
totally dropped as we have discussed form the quarkonium bound state equation. Then the 
BSW wavefunction for quarkonia is reduced to 



,0(2/' ^-l) = exp 



exp 



(4.11) 



which is the exact form as we expected from the qualitative analysis. Here we have not taken 
the limit of mg 00 for heavy quarkonia. Thus a possible mg dependence in wavefunction 
may be hidden in the parameter uj. 

In the following, we shall use the above Gaussian-type wavefunction ansatz to solve the 
light-front quarkonium bound state equation, and from which to determine qualitatively the 
nonperturbative scaling dynamics and to examine the ideas of WCT to heavy hadron bound 
states. 
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C. WCT of nonperturbative QCD description to heavy quarkonia 



We take the normahzed wavefunction ansatz of ( [4. Ill ), 



0oo(l/,«;±) = 4\/a( ^ 



3/4 



exp 



2^1 



exp 
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(4.12) 



as a quarkonium trial wavefunction, where uj\ means that the wavefunction is also scale 
dependent. Substituting the above wavefunction into the quarkonium bound state equation 
( ^!7|) , we have 



— - A L/ /■ 111 c -r >^nonc 



-A Cf In e + Snonc + S< 



Couli 



(4.13) 



where Skin = (Si^a + represents the kinetic energy, Snonc is the contribution of the 
noncancellation of the instantaneous interaction, 



27r2 



CfX^L + In ^ + Ei(tzr2) + y^Eiiizu)], 



(4.14) 



and Scoui from the Coulomb interaction. 



Coul 



9l 



27r2 



IT 



Erf(ti7) 



1 

+ — 



l-e 



(4.15) 



where 7 = 0.57721566... the Euler constant, e the small longitudinal momentum cutoff, the 



dimensionless w is defined by w 



and El and Erf are the exponential integral function 



and the error function, respectively. 

We may rewrite the term In ^ in Enonc by In ^ = In e + In cc^ + In p- . It shows 
that Snonc contains a logarithmic divergence In e which exactly cancels the same divergence 
from the self-energy correction, as expected, and the term In uj"^ is the logarithmic confining 
energy. 

After the cancellation of the infrared Ine divergences in eq.( |4.13| ), the binding energy for 
heavy quarkonia is given by the kinetic energy plus the interaction energy: 



2r 



Skin ~l~ Scon f ~l~ St 



^C'oul 



-I 



3ul + + 



A_ 

A2 



(4.16) 



where 



F{tu) = 7 + InccJ^ + Ei{tu^) - ^ 1 - 2Erf(u7) 



1 






"1 - e-^'' 







(4.17) 



In Fig.l, we plot the confining energy, the Coulomb energy and the totally interaction energy 
as functions of w which is proportional to the radial variable in the light-front space. 



1 



(4.18) 
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Fig.l shows that the total interaction energy is a typical combination of the Coulomb in- 
teraction at short distance and a confining interaction at long distance that has been widely 
used in previous phenomenological description of hadronic states, but it is now analytically 
derived from QCD without introducing any free parameter except for the QCD running cou- 
pling constant. Furthermore, eq.( |4.16| ) also indicates that without considering the kinetic 
energy, we cannot find a stable quarkonium bound state. The kinetic energy balances the 
interaction energy and ensures the existence of a stable solution for ( |4.16| ). Therefore, the 
first order kinetic energy in HQET is an important nonperturbative effect in binding two 
heavy quarks, as noticed first by Mannel et al. |TT]. 



If we know the experimental value of the quarkonium binding energy A, minimizing 
eg. ([4.161) can completely determine the parameter u!\ and the coupling constant g\. The 



precise value of quarkonium binding energy that can be compared with the data in Particle 
Data Group must include the spin-splitting energy [l/mq corrections) which we will 
present in the forthcoming paper p^. Here, to justify whether a WCT of nonperturbative 



QCD can become possible in the present formulation, we will give a schematic calculation. 
It is known that A is of the same order A.qcd which is about 100 ~ 400 MeV. We choose A as 
a typical hadronic mass, A = 1 GeV. The charmed and bottom quark masses used here are 
rric = 1.4 GeV and mh = 4.8 GeV. From the particle data the lowest charmonium states 



M(?7,(15)) = 2.98 GeV and M(J/^(15) = 3.10 GeV^the bottomonium state M(T(1^)) 
9.46 GeV. Hence the binding energies of quarkonia A < 400 MeV in the above choice of 
heavy quark mass. To solve ( [4.16| ) we shall take several values of A within the above range. 
(In principle, it is not necessary to must use a value of A below 400 MeV.) The results are 
listed in Tables I and II for charmonium and bottomonium, respectively, where ujxq denotes 



the minimum point of the binding energy A ( [4.16|) 



We see from the Tables I-II that the coupling constant ax = 9x1 is very small. For 
instance, with A = 200 MeV, we obtain 

_ J 0.02665 charmonium, , . 

~ \ 0.06795 bottomonium, ^ ' ' 

which is much smaller than that extrapolated from the running coupling constant in the 
naive perturbative QCD calculation. The parameter ujxq is the mean value of the (transverse) 
momentum square of heavy quarks inside quarkonia: 

{kl) = (4.20) 

For charmonium, we can see that the resulting ujx are typical values of A.qcd ~ A. The 
kinetic energy is about a half of the interaction energy. For bottomonium, we find that the 



binding energy A cannot be too large. In fact, when A is over about 260 MeV, eq.( 4.16 ) 
has no solution. Meanwhile, compared to charmonium, the effective coupling constant is 
relatively large (in contrast to the perturbative running coupling constant which is smaller 
with increasing mq if it is taken as the energy scale). Also the values of uJx in bottomonium 
wavefunctions are larger than that in charmonium. The difference between charmonium and 
bottomonium in the above nonperturbative calculation can be understood as follows. As we 
know, in the nonrelativistic quark model, the quark momentum in quarkonia is proportional 
to the quark mass, uj\ ~ mq [^. Our relativistic QCD bound state solution exhibits such 
a property. This also indicates the first order kinetic energy has of the same order the 
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leading Hamiltonian while all other terms are suppressed by the power of 1/m [see eq.(3.2)]. 
Therefore, the kinetic energy must be treated nonperturbatively when we apply HQET to 
heavy quarkonia. As a result, the bottomonium kinetic energy (~ uj) becomes large as well. 
To have a nonperturbative balance between the kinetic energy and the interaction energy 
in the bound states, the coupling constant in bottomonium must be larger than that in 
charmonium. All these properties now have been manifested in the solution of eq. ( |4.16| ). A 
more precise determination of ax (i.e., gx) requires an accurate computation of the low-lying 
quarkonium spectroscopy with the 1 / ttlq corrections included |^ . Nevertheless, it has been 



shown that the effective coupling constant in Hx is very small at the hadronic mass scale. 

In order to see how this weak coupling constant varies with the scale A, we take A = 200 
MeV and vary the value of A around 1 GeV. The result is listed in Table III. We find that the 
coupling constant is decreased very faster with increasing A. In other words, with a suitable 
choice of the hadronic mass scale A in SRG, we can make the effective coupling constant ax 
in Hx arbitrarily small, and therefore the WCT of nonperturbative QCD can be achieved in 
terms of Hx such that the corrections from Hxi can be truly computed perturbatively. 



V. DISCUSSION AND SUMMARY 



Thus far, the main ideas of the WCT to nonperturbative QCD proposed in the recent 
publication [|l| have, at least qualitatively, been achieved for heavy quarkonia when A is 
around hadronic mass scale (~ 1 GeV). To have a deep understanding of how WCT works 
to nonperturbative QCD, we shall study the A-dependence of ax and discuss the physical 
implication of the WCT in this last section. 



A. Running effective coupling constant ax in SRG 

\ 2 

From Tables I to III, we find that values of the dimensionless parameter w = — are 
greater than 2.5. When x > 2.5, the exponential integral function and the error function 
are simply reduced to Ei(x) = and Erf (a;) = 1. Then, eq. (|4.16|) is reduced to 



2^ = — ^A-^C';A^ — r^-v^^-(7 + ln--lna;y + — , (5.1) 



mg ^ 27r2 [ A^ A^ V' 4 } mg 

with an error less than 10~^. Minimizing A with respect to Ux, we obtain 

Therefore, eq.(|5.1|) becomes 



where uxo is a solution of i f>.2\ ). Since A is the binding energy of quarkonia, it should be 
A-independent. Thus, eqs. (|5.2|) and (|5.3|) determine the A-dependence of gx- The result is 
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9l - /A^- 

"A = — = 



a' 



A 



2 (for a relatively large A >> A), (5.4) 



where the coefficients a and b can be obtained by numerically solving eqs. (|5.2|) and (5.3). 



The coefficient b is almost a constant (with a weak dependence on tuq but independence on 
A and A), while a depends on both A and mq, and also slightly on A. For A > 0.6 GeV, 
the A-dependence in the parameter a is negligible. In Fig. 2, we plot the A-dependence of 
the effective coupling constant ax for charmonium. The dots are the numerical solutions of 
( |4.16| ) and the solid line is given by the analytical form ( |5.4|) with b = 1.15, and a = —0.25 
for A = 0.2 GeV and a = 1.1 for A = 0.4 GeV. We can see that ( p.4| ) is a very good analytical 
solution of the eqs. (|5.2| ) and (|5.3| ) [or of the minimizing eq. (|4.16|) ]. 



To understand the running behavior of ax, we calculate its SRG /5-function. Denote the 
binding energy 

A = A{gx,uJx,X). (5.5) 

The invariance of the binding energy A under similarity renormalization group transfor- 
mation means that A determined from Hxo and Hy^ must be the same for A 7^ A'. Let 
A' = A + 5A, we obtain the corresponding similarity renormalization group equation 

where the quantity (3 is the similarity renormalization group f3 function which is defined by 



(5.7) 

A=A{gA) 



and 'jcu is an anomalous dimension that describes the running properties of the bound state 
wavefunction. The P function can be determined by solving the SRG equation ( |5?7| ) or 



directly computed from eq. (|5.4| ) 



^ = -^^('+a + 2MnA/A) 



A=A(g;,) 

^ —gx (for a relatively large A >> A). (5.8) 

As we will see later, this /3-function also indicates the existence of confining interaction in 
low energy region. 



B. Why a WCT in low energy QCD region 

As we have pointed out in the last section we can always make the effective coupling 
constant ax small with a suitable choice of the low energy scale A such that a WCT to 
nonperturbative QCD can be realized on the light-front. The above dependence of ax with 
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A shows that the effective couphng constant a\ becomes smaller with larger A. Thus, to 
have a WCT to nonperturbative QCD, A cannot be too small. In practice, A is of the order 
1 GeV. Now, the question is what is the physical picture behind this dependence of ax with 
A. 

As we discussed in the previous sections, the effective low energy Hamiltonian Hxq is 
derived by integrating out the light-front high energy modes above A through the SRG. The 
existence of an explicit quark confining interaction in Hxo (like to extract a mean-field in 
strongly correlated systems) is crucial to have a WCT in low-energy QCD region. Formally, 
as we have seen that the confining interaction in Hxo comes from the noncancellation of 
instantaneous gluon exchange with energy below the scale A. With the larger A, the more 
the instantaneous interaction contributes to Hxo so that the confining interaction becomes 
stronger. A physical interpretation of the low energy scale A, given by Wilson 



IS 



that A corresponds to a up-bounded value of dynamical gluon mass in low-energy region. 
Introducing A in SRG forces the gluon exchange energy in the qq effective interaction of Hxo 
to be never less than A. This is equivalent to say that dynamical gluons in low energy region 
have a mass of order A such that the energy of gluons propagating in qq channel cannot be 
smaller than this dynamical mass. A large A implies a large dynamical gluon mass which, on 
one hand, results in an explicit strong quark confining interaction in Hxo, and on the other 
hand, makes the residual quark-gluon interactions weak. The effective coupling constant ax 
characterizes these residual interactions {Hxi is expansed in terms of ax)- While, the smaller 
A implies to take the smaller dynamical gluon mass in low energy region. Then the explicit 
confining interactions in Hxo is weakened and the residual interactions become stronger (i.e., 
some of the confining interaction is hidden in the residual interactions). The limit A — >■ 
implies to keep massless gluons in the low energy region. Then, our effective QCD theory 
turns back to the canonical form where the confining interaction is completely hidden in the 
beautiful but complicated quark-gluon interactions and the coupling constant becomes very 
strong in low energy region, as we expected from the canonical theory. This is a physical 
interpretation of the above dependence of ax with A. It also provides a physical picture why 
a WCT can be applied to low energy QCD region for a finite A in our description. 

To further examine the above understanding, we study the parameter ujxo (the mean 
value of the transverse momentum quarks carried inside the hadrons) as a function of A. 
The result is plotted in Fig. 3. From Fig. 3, we see that with increasing A, uxo is de- 
creased. Correspondingly, the distance between two quarks inside the quarkonia, r/ ~ 
is increased. This ensures that confining interaction becomes stronger with larger A (as 
also shown in Table IV). Besides, Fig. 2 also shows that A determined from the bound state 
equation can never go to zero. It implies that the dynamical gluons in low energy region are 
most likely massive around a few hundreds MeV. 

It may also be worth noting that the low energy scale A introduced here is totally 
different from the high energy scale used in perturbative QCD. The use of the scale 

in perturbative QCD is to separate high and low energy region such that one should in 
principle integrate out the physics below the energy scale [practically, instead of doing 
this (since none is able to do it), one usually uses factorization theorem to experimentally 
fit low energy physics] and then perturbatively compute the high energy physics above 

from QCD. Here, the procedure seems to be an inverse processes of the perturbative 
QCD treatment. We introduce the low-energy scale A and integrate out the high energy 
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dynamics above A. The left is just the low energy QCD theory which describes the low 
energy physics. Normally, by integrating out high energy modes, the resulting low energy 
QCD is not neccessary to lead to the explicit quark confining interaction we obtained in 
HxQ. Therefore the theory is still nonperturbative. It is only on the light-front that the 
high energy modes contain long distance interactions which come from small longitudinal 
momentum of quarks and gluons [see the light-front dispersion relation k~ = pr(/c^ + m^), 
where a large light-front energy k~ can be obtained from a small k^]. Thus, integrating 
out light-front high energy modes results in an explicit quark confining interactions in Hxq. 
Hence, the light-front formulation of the effective theory derived from SRG may be another 
crucial point for making the WCT to nonperturbative QCD to be possible. 

As we know in the canonical QCD theory, the confining interaction should become more 
important if the scale would be smaller. Correspondingly the running coupling constant 
a(Q^) becomes larger. Compare this fact with the above analysis, there may be an inverse 
correspondence between the effective coupling constant ax in Hx and the running coupling 
constant in the full QCD theory: 

a.~=^, with A^^i,. (5.9) 

In other words, the weak-coupling treatment of the confining Hamiltonian Hx may corre- 
spond to an inverse expansion of the original strong-coupling in the full QCD theory. SRG 
just provides an implicative realization for such an expansion. 

More specifically, the running coupling constant in full QCD theory is given by 

where t = |ln^, and is a space-like momentum (the same as A^). Since the similarity 
renormalization group f3 function of eq.( |5.8| ) is determined in the physical sector of low 
energy QCD dynamics, the low energy f3 function of the running coupling constant 'g{Q^) in 
the full theory should behave qualitatively the same. Then we may take the f3{g) function 
in the above equation the same form as eq.( ^.8| ) for « ^^cd (i^ 1°'^ energy region). 
This leads to 

am = «.(/i')^ = co^, Q' « AlcD, (5.11) 



where cq = as{fi^)fi'^ / Aq^j-,. This is consistent with eq.( |5.4|) via eq. (|5.9| ). To give a qual 



itative determination of the coefficient cq in (|5.11| ), we consider A = 0.75 ~ 1.5 GeV and 
A = 0.2 GeV, then A^ = (14 ~ 56)A^ » A^ From (§3), we have: 

-r-2 

«A = (1.0~1.5) — . (5.12) 

The corresponding ~ << A^ ~ ^qcd- From (]5l9|), it follows that for << Ag^^ 

the running coupling constant in the full QCD theory may be approximately given by 

a(g2) = (1.0~1.5)^. (5.13) 
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This is just a qualitative estimation of the running couphng constant in the full QCD theory 
in low energy region that is obtained by the use of the same /3-function determined in SRG 
from our effective low energy QCD description. 

Up to date, no one precisely knows how the full QCD coupling constant varies in low 
energy region. However, it is interesting to see that the running coupling constant given by 
eq. (|5.11|) for small is indeed the basic assumption of the Richardson QQ potential [^]: 



Vm = -Cf^. (5.14) 



where 
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The Richardson QQ potential is proposed to exhibit the asymptotic freedom of QCD in 
short distance and a confining potential in large distance limit. From eq.( |5.14| ), we see that 
for small Q"^ (Q^ << ^qcd)^ 



Comparing with eqs.( |5.11"D and ( [5. 161 ), we have from the Richardson QQ potential(with 

Nf = 3 m) 
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Co = — r = 1.4 . (5.17) 

" (33 -2A^/) ^ ^ 

This result agrees very well with eq.( |5.13| ). Hence, the dependence of ax with A determined 
in our description is essentially a consequence of exhibiting an explicit quark confining 
interaction in our low energy QCD Hamiltonian. The possible inverse relation between ax 
with as{Q'^) discussed above can be considered as another reason why the present effective 
theory to low energy of QCD can be treated as a weak-couple problem. 

Finally, we should also note that there is an very interesting question remained to be 
answered: where is the cross-point (or a cross region) between low energy and high energy 
domains that let ax and as match at this point (or region). At this point, we have no 
answer for it but this is a very important open question for a complete understanding of 
QCD dynamics 0]. 



C. A possible criterion for the A scale- fixing 

Because of the dependence of ax with A (|5.4| ), the realization of WCT to nonperturbative 
QCD in our description depends on a suitable choice of A, as we have mentioned. One may 
further ask what is a suitable choice of A so that we can always treat low energy QCD as a 
weak-coupling problem, and whether there is any ambiguous in the choice of A. 

Table IV tells us that the confining interaction plays a more important role than the 
Coulomb interaction in the determination of the quarkonium bound states on the light-front. 
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This result is different from the usual understanding in the nonrelativistic phenomenolog- 
ical description that the dominant contribution in heavy quarkonium spectroscopy is the 
Coulomb interaction. This discrepancy can be understood as follows. The currently rela- 
tivistic light-front description for heavy quark systems mostly uses the heavy quark masses 
of rric = 1.3 ~ 1.4 GeV and mi, = 4.7 ~ 4.8 GeV or less. In Particle Data Group 



nib = 1.0 ~ 1.6 GeV and rric = 4.1 ~ 4.5 GeV. Thus, the heavy quarkonium binding energies, 
A = Mh — ^rriQ, might be positive [the lowest charmonium ground state M{ric{lS)) = 2.98 
GeV, and the bottomonium M(T(15')) = 9.46 GeV]. Therefore, the Coulomb energy be- 
comes not important for a relativistic description in quarkonia. The dominant contribution 
for binding quarkonium states must come from the nonperturbative balance between the 
kinetic energy and the confining energy. While, in the nonrelativistic phenomenological 
description, one used larger quark masses, rric > 1-8 GeV and nib > 5.1 GeV such 



that the binding energy is forced to be negative. As a result, the Coulomb interactions 
must be dominant in this picture. Of course, on the light-front, the structure of the bound 
state equation is different from the nonrelativistic Schrodinger equation. There is no direct 
comparison. A real solution to the above discrepancy can be obtained after including the 
spin-splitting interactions (l/mg-corrections). 

Nevertheless, from the above analysis, we have seen that the quark confining energy is 
a dominant contribution to the binding energy of bound states. A large binding energy A 
requires a large contribution from the confining interaction and therefore a large value of 
A is needed. A small A requires a relatively small effect from the confining interaction so 
that a small A can be chosen. In all the cases, the effective coupling constant a\ is kept 
at a small value. Thus, by choosing A to fix the rate y at a certain value, the resulting ax 
can remain to be small. The WCT to low energy QCD dynamics is then always applicable. 
This is a criterion for fixing A in our description, similar to the scale-fixing procedure of 
Brodsky-Lepage-Mackenzie in solving scale ambiguity for the usual perturbation theory 
Of course, in principle one can choose an arbitrary value of A, physics should not be 
changed since the dependence of ax with A is determined by the SRG invariance. The above 
criterion for the scale-fixing is useful for the realization of a WCT to perturbative QCD in 
our description. 

In practice, we also find that there is no much freedom in the choice of A value. The scale 
dependence of wavefunctions also provides a restriction on the range of A. For quarkonia, 
ujx is decreased with increasing A. However, ux is proportional to the mean value of the 
(transverse) momentum square of quarks inside quarkonia, which characterizes the size of 
hadrons. Therefore A should not be too large for the best description of bound states. 

On the other hand, as we argued that gluons are massive at the hadronic scale. Hence, 
the lowest bound value of A can be considered as the constituent gluon mass. As an example, 
one may take A to be a constituent gluon mass (about a half of the glueball masses, such as 
the recent possible evidences of /o(1500) and ^(2230) p2|). 



A : (0.75 ~ 1.5) GeV. (5.18) 
Then the effective couphng constant (with A = 200 MeV) is 

OA = 0.06 -0.01 (5.19) 
which is very small. Only a finite A at the order of hadronic mass scale can effectively 
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turn the nonperturbative contribution in the higher order processes into the long distance 
two-body confining interactions through SRG. For the range of eq.( |5.18 ), we further have, 



cuxo = 0.24 ~ 0.2 GeV, (5.20) 

which leads to 

(r) ~ — = 0.8 ~ 1.0 fm. (5.21) 

This gives a resonable quarkonium size. Thus, we have provided a qualitative criterion and 
a quantitative analysis on the choice of A such that a true WCT to nonperturbative QCD 
can always be realized by H\. 



D. Summary 

In summary, to examine the weak-coupling treatment of nonperturbative QCD recently 
proposed in Ref. [Q, we have studied explicitly the heavy quark bound state problem, based 
on the light-front heavy quark effective theory of QCD Firstly, we have used the 

similarity renormalization group approach |jl|J^ to derive the effective confining Hamiltonian 
in the low energy scale for heavy quarks in heavy mass limit. To make the similarity 
renormalization approach practically manable, we have introduced a local cutoff scheme 
( p. Ill ) to the bare QCD (and the effective heavy quark) Hamiltonian, which simplifies the 
cutoff scheme in [Q. Meanwhile we have also introduced a simple smearing function fxij 
( p.l3| ) to the similarity renormalization group approach which further simplifies the original 
formulation of [^. The resulting low-energy effective QCD Hamiltonian of heavy quark 
interactions exhibits the coexistence of a confining interaction and a Coulomb interaction 
on the light-front without introducing any additional free parameter except for the effective 
QCD coupling constant. 

The weak-coupling treatment can be realized for nonperturbative QCD because the light- 
front similarity renormalization group approach with a finite A extracts an explicit quark 
confining interaction from the higher order nontrivial quark-gluon interactions into Hxq so 
that the residual quark-gluon interactions become weak even in low energy region. The weak- 
coupling treatment of nonperturbative QCD is manable for A being around the hadronic 
mass scale. The similarity renormalization group invariance can remove the A-dependence 
in all the physical observables obtained from the effective Hx- The well-defined bound state 
description in QED is a special case (A 0) of the SRG approach. The whole idea of the 
WCT to nonperturbative QCD is originally motivated from the bound state description of 
QED Now, a consistent connection between QCD and QED and their differences in 
determining bound states is explicitly examined on the light-front. 

The applications of the present theory to heavy quarkonium spectroscopy and various 
heavy quarkonium annihilation and production processes can be simply achieved by numeri- 
cally solving the bound state equations ( [4.7|) , and by further including the l/mg corrections 
(which naturally leads to the spin splitting interactions). The extension of the computa- 
tions to heavy-light quark systems is straightforward. The extension of the present work to 
light-light hadrons requires the understanding of chiral symmetry breaking in QCD which is 
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a new challenge to nonperturbative QCD on the light-front. Nevertheless, the present work 
has provided a preliminary realization to the weak-coupling treatment of nonperturbative 
QCD proposed recently by Wilson et al. The new research along this direction is in 
progress. 
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FIGURES 



FIG. 1. A plot of the confining energy, the Coulomb energy and the total interaction energy as 
functions of the dimensionless variable w but vj is proportional to ~ ri via The energies are 
scaled by the factor -Y^Cf- 



FIG. 2. The A-dependence of the effective coupling constant a\. The solid line is given by 
the analytical result (^.4|), and the dots are obtained by numerically minimizing the quarkonium 



binding energy (4.16). Here A is given in units of GeV. 



FIG. 3. The A-dependence of the wavefunction parameter wao which is the solution of mini- 
mizing the quarkonium binding energy ( P^ ), where A is given in units of GeV. 
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TABLES 



Table I. Solution for charmonium ground state with rrif. = 1.4 GeV 



A (GeV) 


o\ 


c^Ao (GeV) 


£un{G^^) + Spot{GcY^) = 2A'(GeV2) 


0.2 


0.02665 


0.222 


0.026836 


0.053173 


0.080009 


0.3 


0.06480 


0.275 


0.067902 


0.112018 


0.179920 


0.4 


0.11831 


0.314 


0.130225 


0.189781 


0.320006 



Table II. Solution for bottomonium ground state with = 4.8 GeV 



A (GeV) 


Ha 


L^'Ao (Ge\ ) 


^A„,(GeV2) + 8,AGcW^) = 2A'(GeV2) 


0.15 


0.029965 


0.492 


0.023397 


0.021602 


0.044999 


0.20 


0.06795 


0.623 


0.050183 


0.029816 


0.079999 


0.25 


0.1385 


0.779 


0.098074 


0.026930 


0.125004 



Table III. Some numerical solution on the A-dependence 
of the weak coupling constant ax. 





charmonium 


bottomonium 


X (GeV) 


Ha 


L^'Ao (GeV) 




L^'Ao (GeV) 


0.75 


0.05960 


0.241 






1.0 


0.02665 


0.222 


0.06795 


0.623 


1.5 


0.00912 


0.199 


0.01607 


0.478 


2.0 


0.00441 


0.185 


0.00695 


0.427 



Table IV. The A-dependence of various interactions to the binding energy 



A (GeV) 


0.5 


0.75 


1.0 


1.2 


1.4 


1.8 


2.0 


3.0 


Skin (GeV^) 


0.04 


0.031 


0.027 


0.025 


0.023 


0.021 


0.02 


0.018 


Sconf (GeV2) 


0.049 


0.050 


0.053 


0.055 


0.057 


0.059 


0.06 


0.062 


Scoi (GeV2) 


-0.009 


-0.001 


-0.001 


0.0 


0.0 


0.0 


0.0 


0.0 
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Fig.2 





Fig.l 



